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The detection of the theoretically expected dark matter is central to particle physics 
and cosmology. Current fashionable supersymmetric models provide a natural 
dark matter candidate which is the lightest supersymmetric particle (LSP). The 
theoretically obtained event rates are usually very low or even undetectable. So 
the experimentalists would like to exploit the modulation effect. In the present 
paper we study a specific class of non-isothermal models involving flows of caustic 
rings. We find that the modulation effect arising from such models is smaller than 
that predicted by the isothermal models. 
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■ 1 Introduction 

in 

In recent years the consideration of exotic dark matter has become necessary 
in order to close the Universe El. Recent data from the Supernova Cosmology 
Project □ ' □ suggest that the situation can be adequately described by a 
barionic component Qb = 0.1 along with the exotic components tie dm = 0.3 
and fl\ = 0.6 (see also Turner, these proceedings). 

Since this particle is expected to be very massive, m x > 30GeV, and 
extremely non relajtiyistic with average kinetic energy T < lQQKeV, it can be 
directly detected era mainly via the recoiling nucleus. 

Using an effective supersymmetric Lagrangian at the quark level, see e.g. 
Jungman et al El and references therein , a quark model for the nucleon Era and 
nuclear wave functions □ one can obtain the needed detection rates. They are 
. typically very low. So experimentally one would like to exploit the modulation 

of the-jesjent rates due to the earth's revolution around the sun. In our previous 
work El El we found anhanced modulation, if one uses appropriate asymmetric 
velocity distribution. The isolated galaxies are, however, surrounded by cold 
dark matter .jvhich, due to gravity, keeps falling continuously on them from 
all directionsli!]. It is the purpose of our present paper to exploit the results of 
such a scenario. 

2 The Basic Ingredients for LSP Nucleus Scattering 

The differential cross section can be cast in the form!: 

da(u,v) = 2 ^ U by y [(Eg + ^) F 2 {u) + t spin F n (u)} (1) 
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with 



a (^) 2 {A 2 [(f s -f s ^-^) 2 } (2) 
m^v A 



The functions associated with the spin, and the small coherent term, 

Ey, associated with the vector contribution, are not going to be discussed 
further (see oue earlier work El). In the above expression m^r is the proton 
mass, jjb r is the reduced mass, F{u) and Fu{u) are the usual and isovector 
spin nuclear form factors and u = q 2 b 2 /2, with b the harmonic oscillator 
size parameter and q the momentum transfer to the nucleus. The scale is 
set bya = ^(Gftun) 2 — 0.77 x 10~ 38 cm 2 The quantity u is also related 
to the experimentally measurable energy transfer Q via the relations Q = 
Qou, Qo = [ArriNb 2 ]' 1 The needed parameters from SUSY are fg , f s (isoscalar 
and isovector sacalar respectively). The differential detection rate for a particle 
with velocity v and a target with mass m is 

dR= P(0) A m da{u,v)\v\ (3) 
m x Am N y n 1 w 

where da(u, v) is given by Eq. ( One normally assumes p(0) = O.SGeV/cm 3 
as the LSP density in our vicinity. 



3 Convolution of the Event Rate 

In this section we will examine the consequences of the earth's revolution 
around the sun (the effect of its rotation around its axis is expected to be 
negligible) i.e. the modulation effect. 

Following Sikivie we will consider 2 x N caustic rings, (i,n) , i=(+.-) and 
n=l,2,...N (N=20 in the model of Sikivie et al), each of which contributes 
to the local density a fraction p n of the the summed density p of each type 
« = +,-• and has velocity y n = (y nx , y ny , y nz ) , in units of v = 220 Km/s, 
with respect to the galactic center. 

We find it convenient to choose the z-axis in the direction of the motion of 
the the sun, the y-axis is normal to the plane of the galaxy and the x-axis is in 
the radial direction. The needed quantities are takan from the work of Sikivie 
(table 1 of last Ref. E3) by the definitions y n = v n /v Q ,y nz = v n ^/v ,y nx = 
Vnr/vo, y n y = Vnz/ v o ■ This leads to a velocity distribution of the form: 



N 

f(v') = S(v - v y„) (4) 

71=1 
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The velocity of the earth around the sun is given byl3. 

Ve = + v\ = vq + vi(sinax. — cosacos'yy + cosa sirvyi) (5) 

where a is the phase of the earth's orbital motion, a = around second of 
June. In the laboratory frame we havecl v = v — ve 

4 Undirectional Event Rates 

Integrating Eq. (||) we obtain for the total undirectional rate 

R = Rt -7777 [l-h(a, Q mm )cosa] (6) 
P(0) 

where Q m in is the energy transfer cutoff imposed by the detector, and a = 
Also p n = d n /p, p = 53n=i d n (for each flow +,-)■ In the Sikivie 
model EiJ 2p/p(0) = 1.25. In the above expressions R is the rate obtained □ 
by neglecting the folding with the LSP velocity and the momentum transfer 
dependence, i.e. by 

R = P -^^-VW)Vs + ^ + (7) 
m x ArriN c 

and it contains all SUSY parameters except m c hi The modulation is described 
in terms of the parameter h. The effect of folding with LSP velocity on the 
total rate is taken into account via the quantity t. 

We like to stress that it is a common practice to extract the LSP nucleon 
cross section from the the meusured event rates in order to compare with the 
SUSY predictions . In such analyses the factor of t is commonly ommitted. It 
is clear, however, that, in going from the data to the cross section, one should 
divide by t. 

The undirectional differential rate takes the form 

{^)=R^rtT(u)[l-co S aH(u)} (8) 

The factor T(u) takes care of the u-dependence of the unmodulated dif- 
ferential rate. It is defined so that 

duT{u) = 1. (9) 

i.e. it is the relative differential rate. u m i„ is determined by the energy cutoff 

due to the performance of the detector. u max is determined by the escape 

2 

velocity v esc via the relations: u max = max(^t n = 1,2, ...,N. On the other 
hand H(u) gives the energy transfer dependent modulation amplitude (relative 
to the unmodualated one). 
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Figure 1: The quantities T(u) and H(u) entering the differential amplitude. Thick solid line 
corresponds to m c f li = 30 GeV the intermediate thickness line to m c f li = 80 GeV, the fine 
line to m c /!i = 100 GeV. The rest correspond to larger LSP masses and fall on top of each 
other, f 



5 Discussion of the Results and Conclusions 

We have calculated the the total event rates for elastic LSP-nucleus scatter- 
ing including realistic nuclear form factors. We focused our attention on those 
aspects of the problem, which do not depend on the parameters of supersymmc- 
try other than the LSP mass. The parameter R, normally calculated in SUSY 
theories, was not considered in this work. The interested reader is referred to 
the literature i-fop. a reviewu and references therin and, in our notation, to our 
previous workta • EL 

We specialized our results for the target 127 1. We considered the effects 
of the detector energy cutoff, by considering two typical cases Q m in = 10 and 
Qmin = 20 KeV. We assumed that the LSP density in-our vicinity and the 
velocity spectrum is that of caustic rings of Sikivie et alt3. 

The total rates are described in terms of t and h. In TABLE I we show 
how they vary the detector energy cutoff and the LSP mass. 

The parameters T(u) and H(u) entering the differential amplitude are 
shown in Fig. 1. The shape of T(u) is analogous to that of the isothermal 
models except that the maximum occurs at u = 0.0, rather than at u = 
0.1. The function H(u) shows oscillations, which result in a smaller total 
modulation. Another way of understanding how the cancellations arize is to 
note that for some rings y nz > 1, while for others y n z < 1. 

The maximum occurs around the 2nd of December, something already 
noticed by Sikivie et alEj. Furthermore the modulation is small, h — 0.25, i.e. 
a 5% difference betwen the maximum and the minimum (see TABLE I). It is 
a bit smaller than that of the symmetric models, but a lot smaller than that 
predicted by the asymmetric onesQtl, i.e h = 0.46 
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Table 1: The quantities t and h in the case of the target 53 1 for various LSP masses and 
Qmin in KeV (for definitions see text). 











LSP 


mass 


in GeV 






Quantity 


Qmin 


10 


30 


50 


80 


100 


125 


250 


t 


0.0 


1.451 


1.072 


0.751 


0.477 


0.379 


0.303 


0.173 


h 


0.0 


0.022 


0.023 


0.024 


0.025 


0.026 


0.026 


0.026 


t 


10.0 


0.000 


0.226 


0.356 


0.265 


0.224 


0.172 


0.098 


h 


10.0 


0.000 


0.013 


0.023 


0.025 


0.025 


0.026 


0.026 


t 


20.0 


0.000 


0.013 


0.126 


0.139 


0.116 


0.095 


0.054 


h 


20.0 


0.000 


0.005 


0.017 


0.024 


0.025 


0.026 


0.026 
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